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Abstract 

We used our previously implemented GW approximation (GWA) based on the 
all-electron full-potential projector augmented wave (PAW) method to study 
the optical properties of small, medium and large-band-gap semiconductors: 
Si, GaAs, AlAs, InP, Mg2Si, C, and LiCl. The aim being to study the size of 
both local-field (LF) and the quasi-particle (QP) corrections to the calculated 
dielectric function obtained using the local density approximation (LDA) . We 
found that while the QP corrections tend to align the calculated structures 
in the optical spectra with their experimental counterparts, the LF effects 
don't change these peak positions but systematically reduce the intensities of 
the so called Ei and E2 structures in all the optical spectra. The reduction 
of the intensity of the Ei peak worsen the agreement with experiment while 
that of E2 improves it. We then show that the local- field correction improves 
considerably the calculated static dielectric constants of all studied semicon- 
ductors. Because the static dielectric constant is a ground state property, the 
remaining discrepancy with experiment should be attributed to the the LDA 
itself. On the other hand, as expected, the calculation of the static dielectric 
constant using the GW quasiparticle energies and including the LF effects is 
underestimated for all the semiconductors. The excitonic effects should then 
correct for this discrepancy with experiment. 
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I. INTRODUCTION 



lat the GWA is a useful method for calculating the quasiparticle prop- 
^ and many applications using this method are now available. In 



In a preceding paper0 (referred to as I), we have implemented the so called GWaBI 
using the all-electron projected augmented wave and applied it successfully to compute 
the electronic structure of some small, medium, and large-band-gap semiconductors: Si, 
GaAs, AlAs, InP, Mg2Si, C, and LiCl. We have then shown that the GWA accounts for 
most of the discrepancy between the LDA and experiment regarding the energy position of 
the conduction states. The remaining discrepancy with experiment is believed to be due 
either to the way the decoupling of the core and valence electrons is performed or to the 
non-selfconsistent implementation of our all-electron GW method. 

It is then clear t 
erties of materials,!^ 

particular, it has been successfully used to obtain a variety of physical properties, ranging 
form the band-width narrowing in alkali-metals, and their clusters,il'ii to the understand- 
ing of the surface reconstruction of semiconductorsHS, or the orientational disorder and 
the photoemission spectra of solid Ceo-li^ 

Nevertheless, one of the most interesting quasiparticle application is the understanding of 
the many-body effects on the single particle excitation's spectra, such as the optical spectra 
of semiconductors, and it is still relatively less explored. In this regards. The recent "ab- 
initio" inclusion of local-field and excitonic effects in the calculation of some semiconductor 
dielectric functions is a great achievement. 

i-i 

There is, however, no systematic study of 
the trend of the size of these effects with respect to several types of semiconductors. It is 
therefore important to study the trend of the local-field contribution to the intensities of the 
structures observed in the optical spectra of various types of semiconductors. 

To date all local-field calculations are based on pseudopotential methods. In particular, 
Louie, Chehkowsky, and Cohen used the empirical pseudo-potential (EPP) methodic to 
computed the optical properties of Si. Recently Albrecht et al. used the ab-initio PP for 
Sic3, and Gavrilenko and Bechstedt for SiC, Si, and CE3. However, these PP calculations do 
not agree with the EPP§, and do not provide any definite trend for the size of local-field 
contribution as a function of the semiconductor's type. 

In this paper, we use our newly developed all-electron GWA based on the projector- 
augmented wave (PAW) method" to compute the local-field effect on various types of 
semiconductors. The dynamical dielectric function is computed in the random-phase- 
approximation (RPA) with and without local-field effects using both the LDA and quasi- 
particle energies calculated within our GWA method. The aim being a systematic study 
of the local field effects in all these semiconductors in order to understand its contribution 
to the macroscopic dynamical dielectric function, to the electron loss energy spectrum, and 
to the static dielectric constant. The latter being a ground state property, and should, in 
principle, be well described by an LDA calculation including the local-field effects. 

Our paper is organized as follows: In the second section we introduce briefly our method 
of calculation. In the third section we apply it to determine the optical properties of two 
distinct semiconductor groups: some small and medium band gap semiconductors: Si, GaAs, 
AlAs, InP, and Mg2Si, and some large band gap semiconductors or insulators: C and LiCl. 
We then compare our results with available calculations and experiments. We will also 
discuss the computed values of the static dielectric constant using the LDA with and without 
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local-field effects and the corresponding quasiparticle results. This will lead us to discuss the 
excitonic effect and its contribution to the calculated quasiparticle static dielectric constants. 



II. METHOD OF CALCULATION 

A. Quasiparticles within the GW approximation 

In this paper, we extend the projector-augmented wave method(PAW)li for solve the 
Kohn-Sham equations to the determination of the optical properties of semiconductors. To 
make a quantitative comparison to experimental results we correct the LDA eigenvalues using 
the quasiparticle energies. As described in our previous paper§ we can find the excitation 
energies of the system by solving a quasiparticle equation instead of locating the poles of 
the Green's function. The quasiparticle energies En(k) are determined from the electron 
selfenergy operator Il(r, r', £'„(k)): 

(T + Kxt + ^h)^kn(r) + J rfV'S(r, r', E„(k))^kn(r') = E„(k)^kn(r) (1) 

Here, T is the kinetic energy operator, V^xt is the external (ionic) potential, Vh is the Hartree 
potential due to the average Coulomb repulsion of the electrons. In the GWA, E is approx- 
imated by a convolution with respect to the frequency variable of the Green's function, G, 
with the screened interaction W calculated within the RPA. 



B. Dielectric function and local field effect 

1. Inclusion of local field effects at the RPA level 

In a crystal, which possesses lattice translation symmetry, a small electric perturbation 
-Eo(q+ G, u) of wave vector q+ G and frequency u produces responses -E'(q + G', u) of wave 
vectors q + G'. The G and G' being reciprocal lattice vectors. Thus, the dielectric matrix 
describing these responses, is of the form eG',G(ct;^) and it can be written as: 

E{q + G', u;) = Y. e^Ulq. ^)^o(q + G, u) (2) 

G 

An external macroscopic electric field can be viewed as a perturbation of vanishingly small 
wave vector q and, therefore, the screening of the external macroscopic field is given by the 
matrix element eoo(q, t^) of the inverse dielectric matrix. In insulating crystals, this results 
in a formula for the macroscopic dielectric function: 

e{u) = lim — j — — - — (3) 

which can be rewritten as: 

e(u;) = linieo,o(q,c^) - lim eo,G(q, t^)eG!G'(q' ^)eG',o(q, t^) (4) 
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The first term of tliis equation is tlie interband contribution to tlie macroscopic dielectric 
function and tlie second term represents the local field contribution to e. The determi- 
nation of the macroscopic dielectric constant amounts to the computation of the inverse of 
^G,G'{^y tu). Adler and WiseiH have derived, essentially by an extension of the random-phase 
approximation (RPA), an approximation to ec.G' for longitudinal fields, 



( ^-A Stt ^ [/„,i,„q-/^,i,]Ma'"(k,q)[Ma'^(k,q)]^ 

i2|q+ G||q + G'l ^^^^ En{k - q) - ^m(k) + uj + td 



where n and m are the band indices, /„,k is the zero temperature Fermi distribution, Q is 
the crystal volume and Mq" (k,q) are the matrix elements 

Ma™(k, q) = (v&k-qnle-'^^+^^-'-lvl/k^) (6) 

which are calculated as described in our preceding paper I in the context of the GW approxi- 
mation. In this expression, the time dependence of the field was assumed to be e"*"^* and the 
small positively defined constant 6 guarantees that the matrix elements of e{u) are analytic 
functions in the half upper plane. Such a matrix could be separated into an hermitian part 
^g'*g' ^^'^ anti-hermit ian part ie^^^t according to 

CG.G' (q, uj) = e^G' (q^ ^) + ^^g]g' (q> ^) (7) 
with e*-^-* for positive u given by 

egG-(q,^)= E (^la + cIla + Gl ^^'^^^'^^ [M^a%K<DYHu^-{Ecik)-Kik~cD]) (8) 
and e^^-* defined by a Kramers Kronig (KK) transform as 

(1) ^ \ X , 2 p /"^ , /^'eglG'lq,^') /ox 

eG:G'(q. ^) = ^G,G' + -P du ^ — (9) 

vr JO uj — uj 

It should be noted here that the matrix elements of e^"^^ and e*-^-* could be chosen to be real 
if the inversion is contained in the point group of the crystal. The calculation of the head 
element limq^o ^o,oi.^y ^) of the wing elements limq^o clq^ ^) necessitate special care if 
we want to determine the optical properties of semiconductors when the GW approximation 
or the scissors-shift approximation is used to determine the electronic structure. Instead 
of handling numerically limq^Q M^"^ {k, q) / q where the quasiparticle wave functions V'kn 
and the quasiparticle energies En(k) are to be used, it is reasonable to approximate the 
quasiparticle wave function with the LDA wave function, and take the limit analytically^^: 

lmiMo""^(k,q)/g = q.(nk|p|mk)/(e„(k) - e„(k)) (10) 

Here Ink) and e„(k) are the LDA wave functions and energies for band n and wave vector 
k, respectively. Indeed, it was shown by inspection that for Si the LDA and the GW wave 
functions have more than 99% overlap^. If the scissors-shift approximation is used, it is easy 
to show that 
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=(^)^^H=6g)^f^^(..-A) (11) 



-G,G' 



where A defines the rigid shift of the conduction bands with respect to the valence bands. 

In the above formahsm we have neglected the exchange-correlation contribution to the 
dielectric function. The calculation of this contribution amounts basically to the determina- 
tion of the exchange-correlation kernel /Ta^cli", r') = d^Exc/dp{r)dp{r') = d\4c/dp|p(r)5(r— r'), 
where E^c and V^c are the exchange-correlation energy and potential, respectively, and p(r) 
is the charge density at r. The calculation of K^c is much more complicated in all-electron 
than in a pseudo-potential method, since one has to determine the matrix of the Kernel in 
the Fourier space. Indeed, the FFT of the exchange-correlation kernel converges very slowly 
with the number of G-vectors in reciprocal space because of the oscillating nature of the 
charge density in real space. Fortunately, it has been shown that this kernel contribution to 
the total dielectric function is small.lH 



2. numerical details 

The size of the dielectric matrix is critical for the convergence of the optical spectrum. 
We have found that a size of 65 x 65 for all systems studied here is good for the convergence 
of the optical spectra, except for LiCl where the convergence was achieved for a matrix size of 
181 X 181. The imaginary part of each matrix element Cq Q,(q 0, a;) is evaluated in energy 
intervals of 0.1 eV up to 200 eV. Then the 'real part' e*-^-* is deduced via a KK transformation 
defined previously (see Eq. (|^)). The linear tetrahedron methodiH is employed to perform 
the summation over the Brillouin zone which appears in Eq. (|^). We use 8000 k-points 
in the full Brillouin zone to calculate the head element and 1000 k-points to calculate the 
wing elements and the body elements. The hermiticity of e^^ci^ the time reversal 

symmetry and the symmetry properties are used to reduce the number of matrix elements 
to be computed. 



III. RESULTS AND DISCUSSION 

A. Calculated optical spectra with and without local-field effects 

The dynamical dielectric function of all the semiconductors studied here are calculated 
using Eq. (^. The local-field effects are represented by the second term on the right of 
this equation. To test the accuracy of the all-electron PAW method we have computed the 
imaginary-part of the dielectric function of Si and GaAs without local field and compared 
them to the full-potential linear muffin-tin (FPLMTO) resultS. Fig. |^ shows that the 
agreement with the FPLMTO spectra is excellent. These results are interesting because 
they set, for the first time, the standard for an accurate LDA dielectric function of Si 
and GaAs computed by two different all-electron methods. This is encouraging since the 
LMTO method is a state-of-the-art first principles method for electronic structure, and in 
comparison, the PAW formalism is much simpler, but nevertheless the method doesn't loss 
any accuracy. 

Before presenting our calculated optical properties of semiconductors, we would like to 
mention the different ways we have obtained the optical spectra: (1) We have used our 
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LDA calculated band structure to directly compute the optical spectra. (2) We have used 
the so called scissors-operator energy shift to the LDA eigenvalues. The value of the shift 
corresponds to our GWA correction of the LDA direct band gap at the F point. (3) Finally, 
we used the GWA calculated quasiparticle energies across the Brillouin zone. All these three 
types of calculations were produced with and without local-field effects. 

The accuracy of the macroscopic function depends on the convergence of all the elements 
of the microscopic dielectric matrix. As stated in the previous section we have found that a 
matrix of 65 by 65 G-vectors and the use of 200 bands in the interband transitions produce 
a well converged e(a;), except for LiCl where a matrix of 181 x 181 G-vectors is used. Fig. |^ 
shows different elements of the microscopic dielectric function of Silicon versus photon energy 
up to 70 eV. The highest intensity of these elements is at least one order of magnitude smaller 
than the e(ooo),(ooo) element. We compared our results to these of Gavrilenko and BechstedtH 
and found that the agreement with their results is only at the semi- quantitative level. We 
are surprised to find that Gavrilenko and Bechstedt have a relatively large intensity in the 
band gap of their /me(ooo),(iii) and /me(iii),(2oo)- 

Fig. ^ shows the calculated real and imaginary part of the dielectric function of Si versus 
photon energy up to 10 eV with and without LF and QP energy shift. These calculations are 
compared to the experimental results of Aspnes and Studna.a The dashed curve represent 
the difference between the calculated optical spectra using the calculated GW quasiparticle 
energies and the one obtained using the scissors-operator energy shift to the LDA eigenvalues. 
This spectral difference is small justifying the use of the scissors-operator for the calculation 
of the optical spectra of small and medium gap semiconductors. Because of this small change 
of the dielectric function due to the use of the quasiparticle energies, and because of the 
high CPU cost in obtaining the quasiparticle energies across the whole Brillouin zone, all the 
other small and medium gap semiconductors are calculated using only the scissors-operator 
shift. Notice that the agreement concerning the peak positions is fortuitous, because the QP 
energy shift underestimates the direct band gap at F (for details, see paper I). A calculation 
with an energy shift which reproduces the experimental band gap will slightly overestimates 
the peak positions by about 0.3 eV, and this overestimation is valid for all semiconductors 
studied here. It is also worth mentionning that our preliminary calculation of the excitonic 
effects shows that in the case of Si the Ei, and E2 peaks are shifted by about 0.2 eV towards 
lower photon energies, and the agreement with experiment concerning the postions and 
intensities of these peaks is recovered. 

On the other hand, in agreement with the empirical pseudopotential (EPP) calculation of 
Louie, Chelikowsky and Cohen we have found that the LF effects do not change the position 
of the structures present on the optical spectrumil. But the intensities of the so called Ei 
and E2 peaks are reduced, again in agreement with the EPP. Thus the local-field effects seem 
to improve the agreement with experiment regarding the intensity of the E2 peak and the 
structures in the higher energy part, and worsen the agreement with experiment regarding 
the low energy part where the Ei peak is located. It is surprising that our calculations 
do not agree well with the ab-initio PP calculation of Gavrilenko and BechstedlH which 
is supposed to be similar to the EPM calculation. The latter calculation found that while 
the LF underestimates the Ei peak intensity in agreement with our calculation and EPP, it 
overestimates the E2 peak intensity in disagreement with our calculation and with EPP. The 
calculated spectrum of Albrecht et al^, obtained by solving the Bethe-Salpeter equation, 
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using the PP method, in a special hmit where only local-field effects are included, agrees 
only qualitatively with our calculation and the EPP te^ultB. Their E[ structure, which 
is located in energy above the E2 structure, has a large intensity and disagrees with our 
calculation and other ab-initio or EPP calculation^^H. It is then not clear what makes the 
extra reduction of their E[ peak when the excitonic effects are included. 

We have perfomed similar calculations for GaAs, AlAs, InP, Mg2Si, C and LiCl, which 
we present in Figs. ^, ^ ^, |^, ^, and ^ and compare to available experimental resultsSS For 
all these semiconductor the scissors-operator energy shift, corresponding to the difference 
between the quasiparticle and LDA eigenvalues at the P point, is used to produce the 
quasiparticle optical spectra, except for C where we have represented the dielectric function 
obtained using the GW quasiparticle energies at each point in the Brillouin zone. In this 
latter case the difference between the optical spectra, calculated using the GW quasiparticle 
energies and the scissors-operator energy shift, is shown by a dashed curve in Fig. H, and 
it is found to be much larger than in the case of Si. We believe that the same conclusion 
should be valid for LiCl, however due to the cost of the GW calculation and the absence of 
the experimental results we preferred not to perform the calculation with the quasiparticle 
energies. Moreover, the small dielectric constants of these two large gap semiconductors 
indicate that excitonic effects are important. For example, the large discrepancy between 
our calculated spectra of Diamond and experiment can be attributed to these effects. 

In conclusion, the most important trend of the LF on the optical spectra of all types of 
semiconductors studied here, is that: The Eis intensity reduction disagrees with experiment, 
whereas the ii^2's intensity reduction agrees well with experiment whenever available. It is 
interesting to notice that for the insulator LiCl the LF effects seem to reduce substantially 
the intensities of the peaks at the low photon energy. 



B. Electron-energy-loss function 

Figs. 11, 12 show our calculated electron-energy-loss (EEL) functions —Im[t~^{q 



0,cij)]o,o for small-band-gap semiconductors: Si, GaAs, AlAs, InP, and Mg2Si and large- 
band-gap semiconductors: C and LiCl, respectively. The calculation are done within the 
LDA with and without the local- field effect. Whenever possible the calculation is compared 
to available energy-loss spectra. The local-field effects seems to improve the agreement with 
experiment by reducing significantly the intensity of the main peak. The EEL functions of 
C and specially of LiCl are much complicated. The LDA C EEL function has two maxima 
at 31.5 and 34.5 eV and these values are shifted to 31.4 and 35.2 eV, respectively, when the 
LF effects are included. The experimental curve seems to present only one resonance at 32 
eV. This discrepancy, could be easly due to small unaccuracy in the calculated dielectric 
function at these high photon energies. 

We did not calculate —Im[e~^{q = 0,uj)]q q for the quasi-particle energy, because we 
believe that GWA is not valid at high energies and as pointed out in Ref.cJ the plasma 
resonance will be pushed towards higher energies in disagreement with experiment. This 
because the electronic structure at higher energy is most probably much better described 
using the LDA than the GWA because: (1) at these higher energies the scattering of an 
electron with the atomic potential is small. In this respect, these high electronic states can 
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be obtained from an almost free-electron theory. (2) the plasmon-pole model is not valid at 
these high energies. 

Table | shows the values of the maxima of the enepx4oss function compared to the ex- 
perimental results obtained from the EEL experiment aill'E3 and from the measured dielectric 
function^. The free-electron plasma frequency is also shown for comparison. The plasma 
resonance of the EEL spectra of Si and GaAs are in good agreement with experiment and 
with the free-electron plasma frequency. The energy electron-energy-loss function of LiCl is 
too complicated containing many peaks, and in the absence of experimental data we prefered 
not to show the values of these maxima in table |. 

C. The static dielectric constant 

The static dielectric function with or without local-field effects is computed using 
the Kramers-Kroenig relations. The calculations were produced using the RPA dielectric 
function and performing analytically the limit q — 0. Table O presents eoo for all semi- 
conductors studied here and compares them to other calculationsi'ii0&S and to available 
experimental result^!. To illustrate our data and stress the agreement with experiment 
we show in Fig. |T^ all our calculated results versus experiment, except for Mg2Si were we 
are not aware of any available experimental result. A perfect agreement with experiment is 
achieved when the calculated value is on the dashed line. Because the static dielectric func- 
tion is a ground state property, we expect that the calculation with LDA including the LP 
effect should reproduce the experimental results. However, we observe only an improvement 
due to this effect. We conclude that the remaining of the discrepancy with experiment is 
due to the used of the LDA itself instead of the full density functional theory. Notice that 
for large band gap semiconductors, the LDA results including the LP effects are in good 
agreement with experiment. This trend was also observed by other researchers0'0. 

On the other hand, the quasiparticle description of the static dielectric constant, is of 
importance, since it will be directly compared to LDA results. A comparison to experiment 
will show whether the excitonic effect corrections are important or not. In this respect we 
notice that, when the GWA energy shift and LP effects are both included, the QP calculation 
slightly underestimates the static dielectric function for all the semiconductors studied here 
regardless of the size of the band gap or the type of semiconductor. This suggests, effectively, 
the importance of the excitonic effects which are expected to produce a positive contribution 
leading to a better agreement with experiment. This positive contribution arises from the 
reduction of the optical band gap from its GWA counterpart, and to the transfer of the 
force of the oscillator towards lower photon energies. It is however interesting to remark 
that since for large band gap materials (C and LiCl) the static dielectric function within 
LDA including LP effects agrees nicely with experiment, the excitonic contribution should 
cancel out the QP correction. 

IV. CONCLUSION 

We have used our previously implemented GWA within the all-electron projected aug- 
mented wave method (PAW) to study the optical properties of some small, medium and 
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large-band-gap semiconductors: GaAs, AlAs, InP, Mg2Si, C and LiCl. In general, the in- 
clusion of the the quasiparticle (QP) energy shift and the local-field (LF) effects improves 
the agreement with experiment. In particular, the LF effects reduce the intensities of the 
so called Ei and E2 peaks without changing their energy positions. This reduction of the 
peak intensity worsen the agreement for the Ei peak but improves it for the E2 peak. This 
trend is observed for all the studied semiconductors and is found to be in agreement with 
the empirical pseudopotential (EPP) calculation of Louie, Chelikowsky and Cohenil. The 
QP energy shift pushes the calculated peaks towards higher energies in agreement with ex- 
periment. However, because the calculated GWA energy shift does not always produce the 
correct experimental band gaps, the agreement of the peak positions with experiment in the 
case of Si, GaAs, and Mg2Si is fortuitous. A calculation using an energy shift which repro- 
duces the experimental band gap will produce theoretical peaks slightly higher in energy 
compared to experiment. On the other hand, our preliminary calculation of the excitonic 
effects for Si shows that the Ei, and E2 peaks are shifted by about 0.2 eV towards lower 
photon energies, and the agreement with experiment concerning the postions and intensities 
of these peaks is is recovered. Thus, the slight shift of the peaks at higher energy, when the 
experimental band gaps are used, is canceled by the excitonic effects. 

The static dielectric function with or without local-field effects is computed using 
the Kramers-Kroenig relations. The calculation were performed using the RPA dielectric 
function and performing analytically the limit q ^ 0. Because the static dielectric function 
is a ground state property, the calculation using LDA and including the LF effect should 
reproduce the experimental results. In our calculations, for small and medium band gap 
semiconductors, we observe only an improvement due to this effect, and concluded that the 
remaining of the discrepancy is due to the use of the LDA itself instead of the full density 
functional theory. However, for large band gap semiconductors, the LDA results including 
the LF effects are in good agreement with experiment. 

On the other hand, the QP calculation, when both the GWA energy shift and LF effects 
are included, underestimated the static dielectric function for all the semiconductors studied 
here. This suggests the importance of the excitonic effects which are expected to produce a 
positive contribution. 
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TABLES 



TABLE L Influence of LF on the energy position of the plasmon peak of the electron-energy-loss 
spectra. Our calculation is compared to available experimental results and to the free-electron 
plasma frequency (in eV). 



Material LDA LDA-^LF Free electron Expt. 

Si 16^6 16^5 16^6 16.4'^, 16.9'' 

GaAs 16.8 16.4 15.6 14.7^ 

AlAs 16.45 15.8 15.8 

InP 15.5 15.0 14.8 

MgaSi 12.65 12.5 13.0 

C 31.5 and 34.5 31.4 and 35.2 31.2 32^ 



^Ref.§, bRef.§, ^Ref.§ 



TABLE IL Influence of LF and QP shifts on the macroscopic dielectric constant Eoo compared 
to other calculations and to experiment .il 



Material 


LDA 




LDA+LF 




QP shift 


QP shift -FLF 


Expt. 


Si 


13.78 


13.6^13.8'' 


12.39 


12.2^,12.4*^ 


12.04 


10.92 


11.7 


Si 




13.75^, 12.8'= 












GaAs 


14.23 


13.1^ 14.17"^ 


12.78 




11.61 


10.51 


10.9 


AlAs 


10.20 


9.5*^ 


8.93 


8.65^ 


8.55 


7.59 


8.2 


InP 


10.71 




9.55 




8.91 


8.01 


9.6 


Mg2Si 


17.73 




15.22 




15.79 


13.56 




C 


5.94 


5.5<= 


5.54 


5.62^^ 


5.25 


4.94 


5.7 


LiCl 


3.35 




2.84 


2.9'' 


2.82 


2.46 


2.7 


^Ref.i, ^ 


Rem, ^ReiB ''ReiB, 


^RefE, ^Ref.@ 
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FIGURES 




FIG. 1. Calculated LDA imaginary part of the Dielectric function of Si and GaAs versus 
photon energy using PAW method (sohd line) and the FPLMTO method0 (dashed line). The 
agreement between the two calculations is excellent. This sets for the first time the LDA results 
of the dielectric functions of Si and GaAs. 
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FIG. 2. Calculated elements of the real (right column) and imaginary part (left column) of 
the symmetrized microscopic dielectric matrix e(q, lj)g,g' of Silicon for the limit q — > and for 
(G,G') = (000,111), (111,111), and (111,200). 
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FIG. 3. Calculated real and imaginary part of the dielectric function of Si versus photon 
energy. The dot-dashed line is our calculation using a rigid energy shift of 0.6 eV of the LDA 
conduction bands, corresponding to our calculated GW band gap correction at the T point. The 
solid line is the calculation with a rigid energy shift of the conduction bands and including the 
local-field effects. The solid line with open circles is the experimental data@. The dashed curve is 
the difference between the calculation using the quasiparticle energy across the Brillouin zone and 
that using the rigid energy shift. This small difference justifies the use of the scissors-energy shift 
for the calculation of the optical properties. 
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FIG. 4. Calculated real and imaginary part of the Dielectric function of GaAs versus photon 
energy. The dot-dashed line is our calculation using a rigid energy shift of 0.75 eV of the LDA 
conduction bands, corresponding to our calculated GW band-gap correction at the T point. The 
solid line is our calculation using the rigid energy shift and including the local-field effects. The 
solid line with open circles is the experimental spectrumB. 



16 





FIG. 5. Calculated real and imaginary part of the dielectric function of AlAs versus photon 
energy. The dot-dashed line is our calculation using a rigid energy shift of 0.95 eV of the LDA 
conduction bands, corresponding to our calculated GW band gap correction at the T point. The 
solid line is our calculation with a rigid energy shift and including the local-field effects. 
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FIG. 6. Calculated real and imaginary part of the dielectric function of InP versus photon 
energy. The dot-dashed line is our calculation using a rigid energy shift of 0.8 eV of the LDA 
conduction bands, corresponding to our calculated GW band gap correction at the T point. The 
solid line is our calculation with a rigid energy shift and including the local-field effects. 
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FIG. 7. Calculated real and imaginary part of the dielectric function of Mg2Si versus photon 
energy. The dot-dashed line is our calculation using a rigid energy shift of 0.33 eV of the LDA 
conduction bands, corresponding to our calculated GW band gap correction at the T point. The 
solid line is our calculation with a rigid energy shift and including the local-field effects. The solid 
line with open circles is the experimental spectrum0. 
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FIG. 8. Calculated real and imaginary part of the Dielectric function of C versus photon 
energy. The dot-dashed line is our calculation using the quasiparticle energies. The solid line 
is our calculation using the quasiparticle energies and including the local-field effects. The solid 
line with open circles is the experimental spectrumEl. The dashed curve represents the difference 
between the calculations using the quasiparticle energies and the LDA calculation with a GW rigid 
energy shift of 1.9 eV of the conduction states. This rigid energy shift corresponds to the GW 
correction of the band gap at the F point. 
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FIG. 9. Calculated real and imaginary part of the dielectric function of LiCl versus photon 
energy. The dot-dashed line is our calculation using a rigid energy shift of 2.8 eV of the LDA 
conduction bands, corresponding to our calculated GW band gap correction at the F point. The 
solid line is our calculation with a rigid energy shift and including the local-field effects. 
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FIG. 10. Calculated static dielectric function compared to experimental results. The filled 
circles represent the LDA values without local- field effects (LF), the open circles the LDA values 
with LF and the up-triangles the LDA without LF but with an energy shift corresponding to the 
GW correction of the direct band gap at the F point, the empty up-triangles are the LDA values 
with the GW energy shift and the LF (see text). A perfect agreement with experiment is achieved 
when a calculated value is on the dashed line. Notice that when the GW energy shift and the LF are 
included the calculation underestimates the static dielectric constant for all these semiconductors 
regardless of the size of the band gap. This suggest the importance of the excitonic effects which 
are expected to produce a positive correction leading to a better agreement with experiment. 
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FIG. 11. Calculated energy loss function with (solid line) and without local- field effects 
(dashed line) of small and medium band-gap semiconductors: Si, GaAs, AlAs, and Mg2Si compared 
to available experimental result£il (solid line with open circle). 
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FIG. 12. Calculated energy loss function with (solid line) and without local- field effects 
(dashed line) of large band gap semiconductors: C and LiCl compared to available experimental 
result^! (solid line with open circle). 
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